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Leptonic and semileptonic D and Dg decays at B-factories 

L. Widhalm (Belle collaboration) 
Institute of High-Energy Physics, Austrian Academy of Sciences, Vienna 1050, Austria 

Recent measurements of branching fractions, form factors and decay constants of leptonic and 
semileptonic decays of _D(a)-mesons acquired at experiments running at the T(4S') resonance energy 
are reviewed. 
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I. INTRODUCTION 

One of the important goals of particle physics is 
the precise measurement and understanding of the 
Cabibbo-Kobayashi-Maskawa (CKM) Matrix. To in- 
terpret results from B-factory experiments such as 
BaBAR [3| and Belle 0, theoretical calculations of 
form factors and decay constants (usually based on 
lattice gauge theory, see e.g. [1]) are needed. It is nec- 
essary to have accurate measurements in the charm 
sector to check (and allow further tuning of) theoret- 
ical methods and predictions. 

Due to their relative abundance and simplified the- 
oretical treatment, (semi) leptonic decays of D or Dg 
mesons are a favored means of determining the weak 
interaction couplings of quarks within the standard 
model. 

This review concentrates on experimental results for 
such decays achieved at experiments running at the 
T{AS) resonance threshold, namely the BaBAR^ and 
Belle experiments. It uses adapted excerpts from the 
cited Belle and BaBAR publications. 



II. THE EXPERIMENTS 

The BaBAR detector [ij reconstructs charged parti- 
cles by matching hits in the 5-layer double-sided sili- 
con vertex tracker (SVT) with track elements in the 
40- layer drift chamber (DCH), which is filled with a 
gas mixture of helium and isobutane. Slow particles 
which do not leave enough hits in the DCH due to 
the bending in the 1.5-T magnetic field, are recon- 
structed in the SVT. Charged hadron identification is 
performed combining the measurements of the energy 
deposition in the SVT and in the DCH with the in- 
formation from the Cherenkov detector (DIRC). Pho- 
tons are detected and measured in the CsI(Tl) electro- 
magnetic calorimeter (EMC). Electrons are identified 
by the ratio of the track momentum to the associated 
energy deposited in the EMC, the transverse profile 
of the shower, the energy loss in the DCH, and the 
Cherenkov angle in the DIRC. Muons are identified 
in the instrumented flux return, composed of resis- 
tive plate chambers interleaved with layers of steel and 
brass. 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 



(SVD), a 50-layer central drift chamber (CDC), an ar- 
ray of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of 
the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in 
detail elsewhere [2| . Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon 
vertex detector were used for the first sample of 156 
fb~^, while a 1.5 cm beampipe, a 4-layer silicon de- 
tector and a small-cell inner drift chamber were used 
to record the remaining 392 fb~^ 



III. SEMILEPTONIC DECAYS TO SCALAR 
MESONS 

Form factors from D meson semileptonic decay have 
been calculated using lattice QCD techniques ,5, ^,0]- 
In the theoretical description, the differential decay 
width is dominated by the form factor /+((?^) @]i 
where is the invariant mass of the lepton pair. Up 
to order m'j it is given by 



Gp\Vcs(d)\ '^ 
247r3 



where Pk{-k) is the magnitude of the meson 3- 
momentum in the D^ig rest frame. 

In the modified pole model |9| , the form factor /+ is 
described as 



/+(0) 



(l-9V<ic)(l-"p9V"^^ole)' 



(2) 



with the pole masses predicted as m{D*) = 
2.11 GeV/c2 (for D^^^ K+l'u) and m{D*) = 
2.01 GeV/c2 (for n+e-iy). Setting ap = 

leads to the simple pole model 

A model independent description of the form factor 
has been studied in [l3|. The most general expres- 
sions of the form factor /^(g^) are analytic functions 
satisfying the dispersion relation: 



Res{f+) 



dt- 



t-q^ 



(3) 
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The only singularities in the complex t = plane 
originate from the interaction of the charm and the 
strange quarks in vector states. They are a pole, 
situated at the D* mass squared and a cut, along 
the positive real axis, starting at threshold (t+ — 
(m_D + rnxY) for D°K~ production. 

This cut t-plane can be mapped onto the open unit 
disk with center at t = using the variable: 



z{t,to) = 



\/^+ — t — — tp 



(4) ^ 



In this variable, the physical region for the semilep- 
tonic decay corresponds to the small (real) range be- 
tween zbzniax = ±0.051. The z expansion of /+, 



f+{t) oc ^afe(to) z''{t,to), 

k=Q 



(5) 
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is thus expected to converge quickly. 



(GeV^/c^) 



A. D ^ (A'/7r)(e//i)!^(e/^) at Belle 

Belle has measured the absolute branching fractions 
and form factors of K~l^vi and — > 'n~l'^vi 

(l — e,fi) 11], using a novel reconstruction method 
with better resolution than in previous experi- 
ments. The analysis is based on data corresponding 
to a total integrated luminosity of 282 fb~^. 

To achieve good resolution in the neutrino mo- 
mentum and q"^, the D° are tagged by fully recon- 
structing the remainder of the event. The studied 
events are of the type e+e~ — *■ D[l^gD*-~X {Dl~ 
-DgjgTTj}, where X may include additional tt^ , tt", 
or mesons (inclusion of charge-conjugate states 
is implied throughout this paper). The D^J^ is re- 
constructed in the modes Z?*+ — > £)'^7r+, TD+tt^ and 
D*o _> D%°,D°j, with D+/" ^ X-(n7r)++/+ 
{n = 1,2,3}. Each -Dtag and i'tag candidate is sub- 
jected to a mass-constrained vertex fit to improve the 
momentum resolution. The 4-momentum of D*~ is 
found by energy-momentum conservation, assuming a 
D[llD;rX event. Its resolut ion is improved by sub- 

jecting it to a fit of the X tracks and the -D^ag mo- 
mentum, constrained to originate at the run-by-run 
average collision point, while the invariant mass is con- 
strained to the nominal mass of a D*~. Candidates 
for TTj are selected from among the remaining tracks, 
and for each the candidate D2;„ 4-momentum is calcu- 
lated from that of the D*-~ and Tr~ ■ The momentum is 
then adjusted by a kinematic fit constraining the can- 
didate mass to that of the For this fit, the decay 
vertex of the D°^^ has been estimated by extrapolat- 
ing from the collision point in the direction of the D^ig 
momentum assuming the average decay length. 



FIG. 1: Belle: Form factors for (a) D" K't+v, in 
bins of 0.067 GeVVc^ and (b) D° Tf-f+i/, in bins 
of 0.3 GeV^/c^. Overlaid are the predictions of the simple 
pole model using the physical pole mass (dashed), and a 
quenched (yellow) and unquenched (purple) LQCD calcu- 
lation. Each LQCD curve is obtained by fitting a parabola 
to values calculated at specific q^ points. The shaded band 
refiects the theoretical uncertainty and is shown within the 
range of for which calculations are reported. 



Background lying under the D^jg mass peak (i.e. 
fake-Dgj^) is estimated using a wrong sign (WS) sam- 
ple where the tag- and signal-side D candidates have 
the same flavor (I^tag instead of ftag)- A MC study 
(including T(4S') — > BB and continuum {qq, where 
g = c, s, M, d) events [l^, [II]) has found that this 
sample can properly model the shape of background 
except for a small contribution from real D^{g decays 
(« 2%) from interchange between particles used for 
the tag due to particle misidentification. Background 
from fake 5^;^ is subtracted normalizing this shape 
in a sideband region 1.84 — 1.85 GeV/c^, yielding 
56461 ± 309stat ± SaOsyst signal D^^ tags. 

Within this sample of D^i^ tags, the semileptonic 
decay D^^^ K^{Ti^)t~v is reconstructed with 
(tt^ ) and £^ candidates from among the remaining 
tracks. The neutrino candidate 4-momentum is re- 
constructed by energy-momentum conservation, and 
its invariant mass squared, m^, is required to satisfy 
Jto^i < 0.05 GeVVc^. 

Multiple candidates still remain in one third of l)gig 
tags, and in about one quarter of the semileptonic 
sample. In these cases all candidates are saved and 
given equal weights such that each event has a total 
weight of 1. 
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TABLE I: Belle: Yields in data, estimated backgrounds, extracted signal yields and branching fractions, where for the 
latter two, the first uncertainty is statistical and second is systematic; small differences in the numbers are due to 
rounding. 



channel 


full Dgig 








n'^e Us TT 




yield 


95250 


1349 




1333 


152 


141 


fake -Dsig 


38789 


12.6 




12.2 


12.3 


12.5 


semileptonic 


n/a 


6.7 




10.0 


11.7 


12.6 


hadronic 


n/a 


11.9 




62.1 


1.8 


9.7 


signal 


5646f 


1318 




1249 


126 


106 


stat. error 


309 


37 




37 


12 


12 


syst. error 


830 


7 




25 


3 


6 


Branching Fraction (fO ^) 


345 ± 10 ± 19 


345 ± 10 ± 21 


27.9 ±2.7 ± 1.6 23.1 ± 2. e 


i± 1.9 


(e and /i channels. 


average) 


345 ± 7 ± 20 




25.5 ± 1.9 ± 1.6 





The contribution from fake D2;„ in the sample of 

big 

semileptonic decay candidates is estimated using the 
_Dgjg invariant mass WS shape of the D^-^^ tag sam- 
ple, normalized in the previously defined sideband 
region. Backgrounds from semileptonic decays with 
either an incorrectly identified meson or where ad- 
ditional mesons are lost in reconstruction are highly 
suppressed by the good neutrino mass resolution. For 
l)gjg ■K^i.~v the most significant background is 
D° K+t-v amounting to 6% - 8% of the to- 
tal yield. It was estimated using the reconstructed 
K^l~v decays in data, reweighted with the 
(independently measured) probability of kaons to fake 
pions. Smaller backgrounds from — > K*^£~iy and 
D° p^£~i' decays amounting to 0.8% — 0.9% were 
measured by normalizing MC to data in the upper 
sideband region > 0.3 GeV^/c^, which is dom- 
inated by these channels. For D^- K^i~v, de- 

_ sig ' 

cays of K*^t'~v contribute at the level of 

0.5% — 0.8%, measured using a sideband evaluation as 
described above, while background from D'^ tt^I^v 
and — > p^i^u was found to be negligible (< 0.07% 
of the total yield). Background from 13°;^ decays to 
hadrons, where a hadron is mis-identified as a lepton, 
is measured with an opposite sign (OS) sample, where 
the lepton charge is opposite to that of the D*~ slow 
pion. Note that the signal is extracted from the same 
sign (SS) sample. In contrast to the SS sample, the 
OS sample has no signal or semileptonic backgrounds; 
fake Zjgjg are subtracted in the same manner described 
previously. Assigning well identified pion and kaon 
tracks a lepton mass, pure background distribu- 
tions are constructed in both SS and OS, which are 
labelled ff^ and /°^, m = K,n. A fit of the weights 
qk and of the components f^^ and f^^ in the ml 
distribution of the OS data sample is performed, and 
the hadronic background in the SS data sample is cal- 
culated as {aKf^ + ayrf^^), utilizing the fact that the 
hadron misidentification rate does not depend on the 
charge correlation defining SS and OS. The method 



has been validated using MC samples. As the muon 
fake rate is about an order of magnitude larger than 
that for electrons, this background is much more sig- 
nificant for muon modes. The signal yields and esti- 
mated backgrounds are summarized in the upper part 
of Table m 

Efficiencies depend strongly on nx, defined as the 
number of 7r*'^°) and mesons assigned to X (in 
e+e^ I?^*gZ3*jg AT), and are determined with MC; 
differences in the nx distribution between MC and 
data give rise to a further (±1.9 ± 3.9)% correction. 
Applying these corrections, the absolute branching 
fractions (normalized to the total number of D^^^ tags) 
summarized in the lower part of Table |T] are obtained. 

The resolution in of semileptonic decays is found 
to be cr,2 = 0.0145 ± 0.0007stat GeVVc^ in MC sig- 
nal events. This is much smaller than statistically 
reasonable bin widths, which have been chosen as 
0.067 (0.3) GcV^/c^ for kaon (pion) modes, and hence 
no unfolding is necessary. Bias in the measurement of 

that may arise due to events where the lepton and 
meson are interchanged, a double mis-assignment, was 
checked with candidate D^-^^ — > K'^i~iy events and 
found to be negligible. The differential decay width is 
bin-by-bin background subtracted and efficiency cor- 
rected, using the same methods described previously. 

The measured q^ distribution is fitted with 2 free 
parameters to the predicted differential decay width 
dT/dq^ of the pole models with /+(0) being one of 
the parameters, and either rripoie (setting ap — 0) or 
Qfp (assuming the theoretical pole) the other. Bin- 
ning effects are accounted for by averaging the model 
functions within individual q^ bins. The fit to the 
simple pole model yields mpo\c[K~ (i^ v) = 1.82 ± 
0.04,tat ± 0.03,y,t GeV/c2 (xV^df = 34/28) and 
mpoic(7r-£+i/) = 1.97 ± 0.08,tat ± 0.04,y,t GeV/c^ 
(x^/ndf = 6.2/10). While the pole mass for the ttIv 
decay agrees within errors with the predicted value, 
m{D*), the more accurate fit of mpoiciKtv) is sev- 
eral standard deviations below m(D*). In the mod- 
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ified pole model, Up describes this deviation of the 
real poles from the m{D*^^-^) masses. Fixing these 
masses to their known experimental values, a fit of Up 
yields ap{D^ R-f+v) ^ 0.52 ± O.OSgtat ± 0.06syst 
(xVndf = 31/28) and ap{D^ TT'l+iy) = 0.10 ± 
0.21,tat ± O.lOsyst (xVndf = 6.4/10). 

The fitted values for /"'""(O) vary little for the 
different fits, for the modified pole model the re- 
sults are f^{0) = 0.695 ± 0.007stat ± 0.022sy,t and 
/;(0) = 0.624 ± 0.020stat ± 0.030syst. 

The measured form factors f^'^{q'^) are shown in 
Figure [T] with predictions of the simple pole model, 
unquenched @ and quenched 0| LQCD. To obtain 
a continuous curve for /+ from the LQCD values re- 
ported at discrete points, the values were fitted 
by a parabola, which is found to fit well within the 
stated theoretical errors and is not associated with 
any specific model. To quantify the degree of agree- 
ment, a x^/ndf is calculated between this measure- 
ment and the interpolated LQCD curve within the 

range for which LQCD predictions are made. For 
the kaon modes, xV^df is 28/18 (34/23), for the pion 
modes 9.8/5 (3.4/5); correlations induced by the fit 
of the calculated q"^ points to a parabola have been 
considered. 



B. D ^ Keve at BABAR 

This subsection is an adapted excerpt of BABAR's pub- 
lication [3. 

The corresponding BABAR analysis [TJ] is using a 
total integrated luminosity of 75 fb~^ collected dur- 
ing the years 2000-2002. It measures the q^ variation 
and the absolute value of the hadronic form factor at 
= for the decay K~e^i'ei'y)- Normal- 

izing to K~TT^, it also gives a value for its 

branching fraction. In contrast to the Belle analy- 
sis, a semi-inclusive reconstruction technique is used 
to select semileptonic decays with less resolution, but 
much higher efficiency. As a result of this approach, 
events with a photon radiated during the decay 
are included in the signal. 

K'^e^i'ei'j) decays are reconstructed in 
e+e^ cc events where the originates from the 
_D*+ D^TT^. Charged and neutral particles are 
boosted to the center of mass system (cm.) and the 
event thrust axis is determined. The direction of this 
axis is required to be in the interval | cos(^thrust )| < 0.6 
to minimize the loss of particles in regions close to the 
beam axis. A plane perpendicular to the thrust axis 
is used to define two hemispheres, equivalent to the 
two jets produced by quark fragmentation. In each 
hemisphere, pairs of oppositely charged leptons and 
kaons are searched for. For the charged lepton candi- 
dates only electrons or positrons with cm. momen- 



tum greater than 0.5 GeV/c are considered. 

Since the Ve momentum is unmeasured, a kinematic 
fit is performed, constraining the invariant mass of the 
candidate e^K^iy^ system to the D'^ mass. In this fit, 
the D'^ momentum and the neutrino energy are esti- 
mated from the other particles measured in the event. 
The D° direction is taken as the direction opposite to 
the sum of the momenta of all reconstructed particles 
in the event, except for the kaon and the positron as- 
sociated with the signal candidate. The energy of the 
jet is determined from the total cm. energy and from 
the measured masses of the two jets. The neutrino 
energy is estimated as the difference between the to- 
tal energy of the jet and the sum of the energies of all 
reconstructed particles in the hemisphere. A correc- 
tion, which depends on the value of the missing energy 
measured in the opposite jet, is applied to account for 
the presence of missing energy due to particles escap- 
ing detection, even in the absence of a neutrino from 
the D° decay. 

Background events arise from T{AS) decays and 
hadronic events from the continuum. To reduce the 
contribution from BB events, selection criteria ex- 
ploiting the topological differences to events with cc 
fragmentation are used. Background events from the 
continuum arise mainly from charm particles Because 
charm hadrons take a large fraction of the charm 
quark energy, charm decay products have higher aver- 
age energies and different angular distributions (rela- 
tive to the thrust axis or to the D direction) compared 
with other particles in the hemisphere emitted from 
the hadronization of the c and c quarks. Selection 
criteria based on these considerations are applied to 
suppress this kind of background. 

The remaining background from cc-events can be 
divided into peaking (60%) and non-peaking (40%) 
candidates. Peaking events are those background 
events whose distribution is peaked around the sig- 
nal region. These are mainly events with a real D*~^ 
in which the slow 7r+ is included in the candidate 
track combination. Backgrounds from e+e~ annihi- 
lations into light uu, dd, ss quarks and BB events 
are non-peaking. To improve the accuracy of the re- 
constructed D'^ momentum, the nominal D*+ mass 
is added as a constraint in the previous fit and only 
events with a probability higher than 1% are kept, 
resulting in 85260 selected candidates containing 
an estimated number of 11280 background events. 
The non-peaking component comprises 54% of the 
background. Detailed studies have been performed to 
understand corrections (and the connected systemat- 
ics) of various components of the peaking background. 

To obtain the true q^ distribution, the measured 
one has to be corrected for selection efficiency and 
detector resolution effects. This is done using an un- 
folding algorithm based on MC simulation of these ef- 
fects: Using Singular Value Decomposition (SVD) [l3| 
of the resolution matrix, the unfolded q^ distribution 
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TABLE II: BABAR: Fitted values of the parameters corresponding to diflerent parameterizations of f+{q^). The last 
column gives the /NDF of the fit when using the value expected for the parameter. 



Theoretical 


Unit 


Parameters 


X^/NDF Expectations 


ansatz 






[xVndf] 


z expansion 




ai = -2.5 ±0.2 ±0.2 


5.9/7 






aa = 0.6 ±6. ±5. 




Modified pole 




Qpoio = 0.377 ± 0.023 ± 0.029 


6.0/8 


Simple pole 


GeV/c^ 


mpoie = 1.884 ± 0.012 ± 0.015 


7.4/8 2.112 [243/9] 



for signal events, corrected for resolution and accep- 
tance effects, is obtained. This approach provides the 
full covariancc matrix for the bin contents of the un- 
folded distribution. To verify that the (f' variation of 
the selection efficiency is well described by the simu- 
lation, a control sample of K~tt^'k'^ is recon- 
structed as if they were K~ e^Ve events, and indicates 
no significant bias. With a second control sample of 
K~TT~^, the accuracy of the £)° direction and 
missing energy reconstruction for the — > K~e^Ve 
analysis is checked. This information is used in the 
mass-constrained fits and thus influences the recon- 
struction. Once the simulation is tuned to reproduce 
the results obtained on data for these parameters, the 
resolution distributions agree very well. One half of 
the measured variation on the fitted parameters from 
these corrections has been taken as a systematic un- 
certainty. 

Effects from a momentum-dependent difference be- 
tween data and simulated events on the charged lepton 
and on the kaon identification have been found to be 
< 2% and included in the corrections and systematics. 
Care has also be taken to correctly understand radia- 
tive decays where = {po —PkY = {Pe + Pv + P-yY ■ 
Corresponding corrections have been applied and the 
corresponding uncertainties enter in the systematic 
uncertainty evaluation. Toy simulations have been 
used to verify that the statistical precision obtained 
for each binned unfolded value is correct and if biases 
generated by removing information are under control. 

The fit to a model is done by comparing the num- 
ber of events measured in a given bin of q^ with the 
expectation from the exact analytic integration of the 
expression |/+(q^)|^ over the bin range, with 

the overall normalization left free. The summary of 
the fits to the normalized q^ distributions is presented 
in Table [Til As long as the form factor parameters are 
left free in the fit, the fitted distributions agree well 
with the data and it is not possible to reject any of 
the parameterizations. As also observed by Belle and 
other experiments, the simple pole model ansatz, with 
"T-poic = ™D; — 2.112 GeV/c^ does not reproduce the 
measurements. 

Figure [2] shows the product P x $ x /+ as a function 



of z{q^) (defined above), where Px $ is the theoretical 
normalization p^ . which constrains this product to 
unity at z = ^max (equivalent to q^ =0). The data 
are compatible with a linear dependence, which is fully 
consistent with the modified pole ansatz for /+(g^). 




-0.05 -0.025 




-z 



0.025 0.05 



FIG. 2: BABAR: Measured values for P x $ x /+ are plotted 
versus ~z and requiring that P x <I> x /+ = 1 for z = 2max. 
The straight lines represent the result for the modified pole 
ansatz, the fit in the center and the statistical and total 
uncertainty. 



The D° ^ K e'^Ve branching fraction is measured 
relative to the reference decay channel, — > K~-k'^: 



Rd 



PP(i?0 ^ if-e+Z^e)data 



BR{D° K-Tr+) 



data 



(6) 



Using slightly adapted selection criteria, after back- 
ground subtraction there remain 76283 ± 323 events 
of D° K-e+Ve in data. To select P»° -> K~tt+ 
candidates, care has been taken to do this in the 
most similar way possible. After background sub- 
traction and the necessary corrections, there are 



6 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



134537±374 candidates selected in the interval 5{m) e 
[0.142, 0.149] GeV/c2. 

A summary of the systematic uncertainties on Rjj 
is given in [14]. The measured relative decay rate is: 



Rd = 0.9269 ± 0.0072 ± 0.0119. 



(7) 



Using the world average for the branching fraction 
BR{D° K-t:+) = (3.80 ± 0.07)% gives 
BR{D^ -> K-e+i^eh)) = (3.522 ± 0.027 ± 0.045 ± 
0.065)%, where the last quoted uncertainty corre- 
sponds to the accuracy on BR{D^ K~tt^). 

The value of the hadronic form factor at = Q can 
be obtained as 



/+(0) = 




Tool' 



(8) 



where BR is the measured K e'^Vp branching 

fraction, td« = (410.1 ± 1.5) x lO^^^ s [iB] is the £)° 

lifetime and / = S^"-\pK{q^f \f+{q^) I f+{0)\\dq\ 
To account for the variation of the form factor within 
one bin, and in particular to extrapolate the result at 
— 0, the pole mass and the modified pole ansatze 
have been used; the corresponding values obtained for 
/+(0) differ by 0.002. Taking the average between 
these two values and including their difference in the 
systematic uncertainty, this gives 



/+(0) 0.727 ±0.007 ± 0.005 ±0.007, 



(9) 



where the last quoted uncertainty corresponds to the 
accuracy on i?i?(_D" K^n^), t^o and \Vcs\- For the 
z expansion, this corresponds to ao — (2.98 ± 0.01 ± 
0.03 ±0.03) X 10-2. 



IV. SEMILEPTONIC DECAYS TO VECTOR 
MESONS 

The differential semileptonic decay rate of a scalar 
meson to a vector meson, specifically, Df (j/e' ^Vp ,, 
depends on the four variables g^, ^e, 6v and x WA - 
depicted in Fig. [31 

Neglecting the electron mass, the differential decay 
rate as function of these four variables depends in a 
given way on three form factors 



MAq^) 

V{q^) 



^1,2(0) 

1 — /m?^ 
1 — q'^ /my 



(10) 
(11) 



with the pole masses rriA — 2.5 GeV/c^ and my = 
2.1 GeV/c^. Measurements have usually been ex- 
pressed in terms of the ratios of the form factors at 
q^ = 0, namely: 

rv = V{0)/Ai{0) and = ^2(0)^1(0). (12) 
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FIG. 3: Definition of the angles 9e, 9v, and x- 



Based on a prediction by [20|, ry is a constant de- 
pending only on particle masses, 

(mn + rrij,)'^ 
rv = ^—^ = 1.8. 13 



A. ~* (j)eue at BABAR 

This subsection is an adapted excerpt of BABAR's pub- 
lication [1^ . 

BaBar has presented a study of the hadronic form 
factors for the vector meson decay Df <f>e^i^e with 
(j) K^K~ [1^ (results still preliminary). This anal- 
ysis is based on a fraction of the total available BaBAR 
data sample, corresponding to integrated luminosities 
of 78.5 fb"^ recorded on the T{AS) resonance. It fo- 
cuses on semileptonic decays of D'^ mesons which are 
produced via e'^e^ cc annihilation. Dg mesons pro- 
duced in BB events are not included and treated as 
background. 

Similar to BaBAR's I?" — + Kei'e analysis presented 
above, a plane perpendicular to the thrust axis is 
used to define two hemispheres, equivalent to the two 
jets produced by quark fragmentation. In each hemi- 
sphere, decay products of the Df, a charged lepton 
and two oppositely charged kaons are searched for. 
Charged leptons are required to have a cm. momen- 
tum larger than 0.5 GeV/c. The unmeasured neu- 
trino momentum is determined in a way similar to 
the D'^ Kevp analysis presented above. 

Figure m shows the K'^K' invariant mass distribu- 
tion for the selected decays compared to MC simula- 
tion and the composition of the background. (f> can- 
didates are defined as K~^K~ pairs with an invari- 
ant mass in the interval from 1.01 and 1.03 GeV/c^. 
Various selection criteria are applied to suppress the 
background About 71% of the total background 
include a true (j) decay combined with an electron from 
another source, namely B meson decays (41%), charm 
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FIG. 4: BABAR: K K~ invariant mass distribution from 
data and simulated events. MC events have been nor- 
malized to the data luminosity according to the different 
cross sections. The excess of signal events in the 4> region 
can be attributed to a different production rate and decay 
branching fraction of Z)J mesons in data and in simulated 
events. Dedicated studies have been done to evaluate the 
amount of peaking background in real events. 



particle decays (25%), photon conversions or Dalitz 
decays (24%), and the rest are fake electrons. These 
(j) mesons are expected to originate from the primary 
vertex, or from a secondary charm decay vertex. 

A maximum likelihood fit is performed to the four- 
dimensional decay distribution in the reconstructed 
variables q^, cos(0y)r, cos{6e)r and Xr using the hke- 
lihood function 



C 



(14) 



i=l 



In this expression, for each bin i, V(ri;. 



datal^MC 



^) is the 

Poisson probability to observe events, when nf^^ 
are expected. Considering the typical resolutions and 
the available statistics, the four variables are divided 
into 4 bins each, corresponding to a four-dimensional 
array with a total number of bins of 625. 

To determine the expected number of signal events, 
a dedicated sample of signal events is generated in 
MC with a uniform decay phase space distribution, 
and each event is weighted using the differential decay 
rate divided by p^. Two of the four variables, cos(6'y) 
and X, are integrated (averaged), taking advantage of 
the fact that the estimated background rate is flat 
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FIG. 5: BABAR: Results from previous experiments and 
present measurement of r2 and rv in Df — > (j>e^i^e decays. 
The error bars represent the statistical and systematic un- 
certainties added in quadrature. These measurements for 
Ds 4>e^Ve decays are compared with the average of sim- 
ilar measurements obtained for D — > K*e'^Ve decays. The 
± one sigma range is indicated by the two parallel lines. 



in these variables. The background components are 
normalized to correspond to the expected rates for 
the integrated luminosity of the data sample. The 
absolute normalization for signal events {Ns) is left 
free to vary in the fit. In each bin («), the expected 
number of events is evaluated to be: 



(0 

m 

(0 

m 



BaBar 
preliminary 



bckff. 



l^tot(Afc) 



(15) 



Here n^'^"**' refers to the number of simulated signal 
events, with reconstructed values of the four variables 
corresponding to bin i. The weight Wj is evaluated 
for each event, using the generated values of the kine- 
matic variables, thus accounting for resolution effects. 



N" 



, (Afe) is the sum of the weights 



for all simulated signal events which have been gen- 
erated according to a uniform phase space distribu- 
tion. Ns and Afc are the parameters to be fitted. 
Specifically, the free parameters A^ are ry, r2, and 
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parameters which define dependence of the form 
factors. To avoid having to introduce finite ranges for 
the fit to the pole masses, m^, we define = 1 + A^. 
This expression ensures that nii is always larger than 
g^_c.0.9GeV^ 

The fit to the four-dimensional data distribution is 
performed using simulated signal events generated ac- 
cording to a uniform phase space distribution. Sig- 
nal MC events are weighted to correct for differences 
in the quark fragmentation process between data and 
simulated events. 

Differences between data and MC have been mea- 
sured using —^ (j)TT~^ decays, according corrections 
were applied. The influence of combinatorial back- 
ground has been studied and considered in the sys- 
tematics. Background from (j) mesons produced in D 
or B decays results in a further correction in MC (to 
calibrate the production rate) and corresponding 
systematics. The effect of uncertainties due to finite 
MC statistics and background estimation on the fit 
has been studied with toy simulations, and was also 
included in the final results. Remaining detector ef- 
fects have been determined with control data samples. 
Using D*+ D°7r+ and D° K-tt+tt° events it 
has been verified that differences between data and 
simulated events in the resolution of the variables 
and cos(0e) are small compared with other sources of 
systematic uncertainties. They have been neglected 
at present. 

Using fixed values for the pole masses (ruA = 
2.5 GeV/c2 and my = 2.1 GeV/c^), the final fit re- 
sults including all corrections are: 

Ns ^ 12886 ± 129 
rv = 1.636 ±0.067 ±0.038 
ra 0.705 ± 0.056 ±0.029. 

Keeping my fixed, for which there is no sensitivity, 
and adding uia as additional free parameter, the fit 
results in 

Ns = 12887 ± 129 
rv = 1.633 ±0.081 ±0.068 
r2 = 0.711 ±0.111 ±0.096 

niA = ^.bSto lt ± 0-54 GeV/c^ 

The measurements of the parameters ry and r2 for 
the semileptonic decay —^ (jye'^i'e have an accuracy 
similar to the one obtained for D K*e^Ve decays 
see Fig. [3 

V. LEPTONIC DECAYS 

The purely leptonic decay 13+ i^i^e (the charge 
conjugate mode is implied throughout this paper) is 
theoretically a rather clean decay; in the Standard 



Model (SM), the decay is mediated by a single virtual 
W^+-boson. The decay rate is given by 

TiDt - £+.,) = ^4,m|Af^. (^1 - 

(16) 

where Gp is the Fermi coupling constant, mi and 
Mjj+ are the masses of the lepton and of the Dg 
meson, respectively. Vcs is the corresponding CKM- 
matrix element, while all effects of strong interaction 
are accounted for in the decay constant /^+ . Due 
to helicity conservation, the decay rate is highly sup- 
pressed for electrons. Since the detection of r's involve 
additional neutrinos, the muon mode is experimen- 
tally the cleanest and most accessible mode. 



A. Ds ^ niy^ at BABAR 

This subsection is an adapted excerpt of BABAR's pub- 
lication [2^ . 

BaRAR performed a measurement of the ratio 
r{Ds fjLVfj) /T{Ds (f>TT) and the decay constant 
, based on a total integrated luminosity of 230.2 

In order to measure Df —f fi^iy^, the decay chain 
D*+ jDf,Df ^J-^i^fi is reconstructed from £)*+ 
mesons produced in the hard fragmentation of con- 
tinimm cc events. The subsequent decay results in 
a photon, a high-momentum and daughter muon 
and neutrino, lying mostly in the same hemisphere 
of the event. Signal candidates are required to lie in 
the recoil of a fully reconstructed D^, D+, Df, or 
D*~^ meson (the "tag") reconstructed in a variety of 
modes [l^ wherein the tag flavor is uniquely deter- 
mined. To eliminate signal from B decays, the mini- 
mum tag momentum is chosen to be close to the kine- 
matic limit for charm mesons arising from B decays. 

For each event a single tag candidate is chosen and 
then used in the subsequent analysis. To pick this tag 
among multiple candidates within an event (there are 
1.2 candidates on average in events with at least one 
candidate) modes of higher purity are preferred. In 
events where two tag candidates are reconstructed in 
the same mode, the quality of the vertex fit of the 
D meson is used as a secondary criterion. After sub- 
tracting combinatorial background there are 5 * 10^ 
charm tagged events with a muon amongst the recoil- 
ing particles. 

The signature of the decay D*~^ iDf is a nar- 
row peak in the distribution of the mass difference 
AM = Mifivj) - Mifiv) at 143.5 MeV/c^ The D*+ 
signal is reconstructed from a muon and a photon can- 
didate in the recoil of the tag. Muons are identified as 
non-showering tracks penetrating the IFR. The muon 
must have a momentum of at least 1.2 GeV/c in the 
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center-of-mass (CM) frame and have a charge consis- 
tent with the tag flavor. Clusters of energy in the 
EMC not associated with charged tracks and exceed- 
ing an CM energy of 0.115 GeV are identified as pho- 
ton candidates. 

The CM missing energy (E^^^^) and momentum 
(Pmiss) calculated from the four- momenta of the 
incoming e~^e~ , the tag four-momentum, and the four- 
momenta of all remaining tracks and photons in the 
event. The energy of the charged particles that do not 
belong to the tag is calculated from the track momen- 
tum under a pion mass hypothesis. Assigning a mass 
according to the most likely particle hypothesis has 
negligible effect on the missing energy resolution. 

The neutrino CM four- momentum (p* — 
is estimated from the niuon CM four-momentum (p* ) 
and p^iss, using a technique adopted from Ref. [23[. 
The difference Ip^jigg — | is minimized, while the in- 
variant mass of the neutrino-muon pair is required to 
be the known mass of the Df [23] ■ The muon CM 
four-momentum (p* ) is combined with p* to form the 
candidate. The candidate is then combined 
with a photon candidate to form the D*^. The se- 
lection requirements on E*^-^^^, Pcom and other vari- 
ables have been optimized to maximize the signifi- 
cance s/Vs + b, where s and b are the signal and back- 
ground yields expected in the data set. 

One class of background are events e+e^ ff, 
where / = u,d, s, b, or r, which do not contain a 
real charm tag. The contribution of these events is 
estimated from data using the tag sidebands. In ad- 
dition there are events e"^e~ — > cc where the tag is 
incorrectly reconstructed. Although these events po- 
tentially contain the signal decay, they are also sub- 
tracted using the tag sidebands. These two sources 
amount to « 42 % of the background. A second class 
of background events (« 26 %) are correctly tagged cc 
events with the recoil muon coming from a semilep- 
tonic charm decay or from r+ — s- ^^v^Vt. This in- 
cludes events jDf — > 7r+i^^, t+ ^i^v^Dr- 
To estimate the size and shape of this background 
contribution, the analysis is repeated, substituting 
a well-identified electron for the muon. Except for 
a small phase-space correction, the widths of weak 
charm decays into muons and electrons are assumed 
to be equal. QED effects such as bremsstrahlung 
(e"*" — > 76^) energy losses and photon conversion 
(7 e+e~), where the muon equivalents have a 
much lower rate, are explicitly removed. In partic- 
ular, bremsstrahlung photons found in the vicinity 
of an electron track are combined with the track. 
The small number of events with an electron from a 
converted photon that survive the selection are sup- 
pressed by a photon conversion veto, using the vertex 
and the known radial distribution of the material in 
the detector. The muon selection efficiency as a func- 
tion of momentum and direction is measured using 
e^e^ — )■ /i+/i~7 events, while radiative Bhabha events 
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FIG. 6: BABAR: AM distribution of charm-tagged events 
passing the signal selection. The tag can be from the tag 
signal region (solid lines) or the sidebands (dashed lines). 
In the bottom plot the signal muon is replaced with an 
electron to estimate the semileptonic charm and r decay 
background. 



are used to quantify the electron efficiency. The ratio 
of muon to electron efficiencies is applied as a weight 
to each electron event. The remaining backgrounds 
are estimated from simulation. 

Events that pass the signal selection are grouped 
into four sets, depending on whether the tag lies in the 
signal region or the sideband regions, and on whether 
the lepton is a muon or an electron (Fig. [6|). For 
each lepton type the sideband AM distribution is sub- 
tracted. The electron distribution, scaled by the rela- 
tive phase-space factor (0.97) appropriate to semilep- 
tonic charm meson decays and leptonic r decays is 
then subtracted from the muon distribution. The 
resulting AM distribution is fitted with a function 
(A^Sig/sig + A'^Bkgd/Bkgd)(AM), where /sig and /ekgd 
describe the simulated signal and background AM 
distributions. The function /gig is a double Gaussian 
distribution. The function /ekgd consists of a double 
and a single Gaussian distribution describing the two 
peaking background components, and a function (2^ 
describing the flat background component. The rela- 
tive sizes of the background components, along with 
all parameters except A'^gig and A^Bkgd are fixed to the 
values estimated from simulation. The fit yields 
A'^sig = 489±55(siai) signal events and has a fit prob- 
ability of 8.9% (Fig.©. 



The branching fraction of 



not be determined directly, since the production 
rate of dI*''^ mesons in cc fragmentation is un- 
known. Instead the partial width ratio T{Df 
/i+t'p)/r(I?+ (j)Tr~^) is measured by reconstructing 
£>*+ ^ j(j)TT+ decays. The D+ fi+iy,, 

branching fraction is evaluated using the measured 
branching fraction for — > (f)Tr~^ . 

Candidate (j) mesons are reconstructed from two 
kaons of opposite charge. The 4> candidates are com- 
bined with charged pions to form Df meson can- 
didates. Both times a geometrically constrained fit 
is employed, and a minimum requirement on the 
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FIG. 7: BABAR: AM distribution after the tag sidebands 
and the electron sample are subtracted. The solid line is 
the fitted signal and background distribution (A^sig/sig + 
A^skgd/skgd), the dashed line is the background distribu- 
tion (A^Bkgd/Bkgd) alone. 
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FIG. 8: BABAR: AM distribution of selected D*+ 
jDf 7(/)7r^ events after the tag sideband is subtracted. 
The solid line is the fitted signal and background distri- 
bution {Ncj,Trf4,Tr + Af07rBkgd/,^.7rBkgd), the dashed line is the 
background distribution (A^,^.n-Bkgd/<i>7rBkgd) alone. 




fit quality is made. The and the Z?+ candidate 
masses must lie within 2 cr of their nominal values, 
obtained from fits to simulated events and data. Pho- 
ton candidates are then combined with the to 
form D*'^ candidates. The same requirements on 
the CM photon energy and D*'^ momentum as in 
the — > jJi^v^ signal selection are made. Data 
events that pass the selection are grouped into two 
sets: the tag signal and sideband regions. After the 
tag sideband has been subtracted from the tag sig- 
nal AA/ distribution, the remaining distribution is 
fitted with (iV^Tr/^TT + A^07rBkgd/07i-Bkgd)(AM), where 
/^TT is a triple Gaussian, describing the simulated 
— > 7D+ 707r+ signal, and /^TrBkgd consists of 
a broad Gaussian centered at 70 MeV/c^ and a func- 
tion [l^l describing the simulated background AA/ 
distributions. The Gaussian describes the background 
T^'^Df — > 7r°(/)7r+ where the photon candidate 
originates from the tt". The relative sizes of the back- 
ground components, along with all parameters except 
N^TT, -^0irBkgd, and the mean of the peak are fixed 
to the values estimated from simulation. The fit 
yields TV^^ = 2093 ± 99 events and has a probability 
of 25.0% (Fig. il). From simulation 48 ± 23 events 
D*+ jDf jfo{980){K+K-)TT+ are expected 
to contribute to the signal, where the error is mostly 
from the uncertainty in the £>+ fo{980){K^K^)Tr^ 
braching ratio. 

Precise knowledge of the efficiency of reconstruct- 
ing the tag is not important, since it mostly cancels 
in the calculation of the partial width ratio. How- 
ever, the presence of two charged kaons in Df — > (/)tt^ 
events leads to an increased number of random tag 
candidates, compared to Df t'-'^'^fj. events, which 
decreases the chances that the correct tag is picked. 
The size of the correction for this effect to the effi- 
ciency ratio {e^^/esig) is determined to be —1.4% in 
simulated events. 

To measure the effect of a difference between the 
D*^ momentum spectrum in simulated and data 



events, D*~^ iDf — + 707r+ events are selected in 
data with the D*~^ momentum requirement removed. 
The sample is purified by requiring the CM momen- 
tum of the charged pion to be at least 0.8 GeV/c. 
The efficiency-corrected D*'^ momentum distribution 
in data is compared to that of D*^ in simulated 
Z?*+ — > "fDf — !■ "f(f)TT^ events. A harder momentum 
spectrum is observed in data. The detection efficien- 
cies for signal and D*"*" iDf ^(pir^ events are re- 
evaluated after weighting simulated events to match 
the D*'^ momentum distribution measured in data. 
The correction to the efficiency ratio is -1-1.5 %. 

With both corrections applied, the partial 
width ratio is determined to be r^^/r07r = 
(iV/e)sig/(Af/e)0. x 6(0 ^ K+R-) 
0.143 ± 0.018(stat), with 6((/) K+R-) = 49.1% 
11. 

A detailed discussion of the systematics can be 
found in ^23| . Using the BaBAR average for the branch- 
ing ratio 6(i:'+ ^ 07r+) = (4.71 ± 0.46) % [H [11], 
the branching fraction B{Df — > l^^v^) — (6.74 ± 
0.83 ± 0.26 ± 0.66) x 10"^ and the decay constant 
fo, = (283 ± 17 ± 7 ± 14) MeV are obtained. The 
first and second errors are statistical and system- 
atic, respectively; the third is the uncertainty from 
B{Dt ^ (/.7r+). 



B. Da fiffj, at Belle 

The Belle analysis uses data corresponding to 548 
fb^^ to study the decay M^^M' using the full- 

reconstruction recoil method first established in the 
study of semileptonic D mesons described above. 

This analysis uses fully reconstructed events of the 
type e+e' D*D'^'°R'^'°X , where X can be any 
number of additional pions from fragmentation, and 
up to one photon _31]. The tag side consists of a D- 
and a R meson (in any charge combination) while the 
signal side is a D* meson decaying to Dgj. Recon- 
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structing the tag side, and allowing any possible set 
of particles in X, the signal side is reconstructed in 
the recoil, using the known beam momentum. 

Tracks are detected with the CDC and the SVD. 
They are required to have at least one associated hit 
in the SVD and an impact parameter with respect 
to the interaction point in the radial direction of less 
than 2 cm and in the beam direction of less than 
4 cm. Tracks are also required to have momenta in 
the laboratory frame greater than 92 MeV/c. A like- 
lihood ratio for a given track to be a kaon or pion, 
which is required to be larger than 50%, is obtained 
by utilising specific ionisation energy loss measure- 
ments made with the CDC, light yield measurements 
from the ACC, and time-of-flight information from the 
TOF. Lepton candidates are required to have momen- 
tum in the lab frame larger than 500 MeV/c. For 
electron identification we use position, cluster energy, 
shower shape in the ECL, combined with track mo- 
mentum and dE/dx measurements in the CDC and 
hits in the ACC. For muon identification, we extrap- 
olate the CDC track to the KLM and compare the 
measured range and transverse deviation in the KLM 
with the expected values. Neutral pion candidates are 
reconstructed using photon pairs with invariant mass 
within ±10 MeV/c^ of the nominal 7r° mass. Neutral 
kaon candidates are reconstructed using charged pion 
pairs within ±30 MeV/c^ of the nominal if ° mass. 

Tag-side D mesons (both charged and neutral) are 
reconstructed in Z? — s- Km: with n = 1,2,3 (total 
branching fraction w 25%). Mass windows have been 
optimised for each channel separately, and a com- 
bined mass- and vertex fit (requiring a confidence level 
greater than 0.1%) is applied to the D meson to im- 
prove the momentum resolution. DJ-candidates are 
not directly reconstructed, but searched for in the re- 
coil of DKX , using the known beam momentum, by 
applying a mass window cut of ±150 MeV/c^ around 
the nominal D* mass ^16] . Since at this point in the re- 
construction X can be any set of remaining pions and 
photons, there are usually a large number of combina- 
torial possibilities. This number is reduced by requir- 
ing the presence of a photon that is consistent with 
the decay £)* Dgj where the Dg has its nominal 
mass within a mass window of ±150 MeV/c^. Further 
selection criteria are applied on the momentum of the 
primary K meson in the e'^e^ rest frame, pK, which 
should be smaller than 2 GeV/c; the momentum of 
the D meson in the e+e~ rest frame, po, should be 
larger than 2 GeV/c; the momentum of the D* me- 
son in the e+e^ rest frame, pd*, is required to be 
larger than 3 GeV/c and the energy of the photon in 
D* — > I?s7, E^,Ds ill the lab frame, is required to be 
larger than 150 MeV/c^, irrespective of 9^. To further 
improve the recoil momentum resolution, inverse [33] 
mass-constrained vertex fits are then performed for 
the D* and Dg, requiring a confidence level greater 
than 1%. After all these selections are applied, the 



average number of combinatorial reconstruction pos- 
sibilities is w 2 per event. The sample is further di- 
vided into a right- (RS) and wrong-sign (WS) part, 
depending on the relative charges of the primary K 
meson, the D meson and that of the K meson the D 
meson decays into {Ko), compared to the charge of 
the D* meson, which is fixed by the total charge of the 
X assuming overall charge conservation for the event. 

Within this sample of Ds-tags, decays of the type 
Ds fJ.Vfi are selected by requiring another charged 
track that is identified as a muon and has the same 
charge as the Dg candidate. No additional charged 
particles are allowed in the event, and additional en- 
ergy corresponding to neutral particles is required to 
be smaller than l/n GeV where n is the number of 
additional neutral particles. After these selections, in 
almost all cases only one combinatorial reconstruction 
possibility remains. Figure [5] shows the mass spectra 
for D^-tags and neutrino candidates. 

nx is defined as the number of primary particles 
in the event, where primary means that the particle 
is not a daughter of any particle reconstructed in the 
event. The minimal value for nx is three correspond- 
ing to a e+e" D*DK event without any further 
particles from fragmentation. The upper limit for nx 
is determined by the reconstruction efhciency; Monte- 
Carlo (MC) shows that the number of reconstructed 
signal events is negligible for nx > 10. As the ef- 
ficiency very sensitively depends on nx, it is crucial 
to use MC that correctly reflects the nx distribution 
observed in data. Unfortunately, the details of frag- 
mentation processes are not very well understood, and 
standard MC |12i| shows notable differences compared 
to data. Furthermore, the true (generated) value 
differs from the reconstructed n^, as particles can be 
lost or wrongly assigned. Thus the measured (recon- 
structed) nx distribution has to be deconvoluted so 
that the analysis can be done in bins of to avoid 
bias in the results. 

To extract the number of Z^^-tags as function of 
in data from the background, 2-dimensional his- 
tograms in n^ (ranging from 3 to 8) and the invariant 
recoil mass tud^ are used. The signal shapes for differ- 
ent values of (ranging from 3 to 9 ISJ) o f the signal 
are modelled with generic MC (CMC) [ij], which has 
been filtered on the generator-level for events of the 
type e+e~ D*DKX. The weights of these com- 
ponents, w^' ,i ^ 1..6, are free parameters in the fit 
to data. As a model for the background in RS, the 
WS data sample is used. Each slice of n^ is fitted 
separately, adding another 6 free parameters. Since 
the WS-sample contains some signal (« 10% of the 
RS signal), these signal components (in shces of n^) 
are also included in the fit as independent parame- 
ters (yielding a negative weight to compensate for the 
WS signal present in the data shapes). The fit is per- 
formed simultaneously with all these free parameters. 
As a crosscheck, the fit has also been performed us- 
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The red shaded areas show the fitted background, the cyan shaded bands show the fit with systematic uncertainties. 
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ing generic MC RS-sample backgrounds, which gives 
a negligible change in the results. A further cross- 
check involved dividing the MC sample randomly into 
two halves, using the shapes of the first half to fit the 
signal in the second. The result as function of ri^, 
normalized to the amount of signal in the first half, 
fits to a flat line as 0.990 ± 0.046, which agrees well 
with the expectation of 1. The total number of recon- 
structed Z?s-tags in data is calculated as 

iVl?>E-^<^^'\ (17) 

where tv^'^*""'* represents the total number of recon- 
structed filtered GMC events that were generated in 
the i-th bin of (regardless of the reconstructed n^ ) 
and wf" is the fitted weight of this component. 

To fit the number of Dg — + /ii/-events as function of 
n^, 2-dimensional histograms in and the missing 
mass squared rn^ are used. The shape of the signal 
is modelled with signal MC. As MC studies show, the 
background under the /xi/-signal peak consists primar- 
ily of non-Dg decays, semileptonic decays (where 
the additional hadrons have low momenta and remain 
undetected) and leptonic r decays (where the r decays 
to a muon and two neutrinos). Hadronic Dg decays 
(with one hadron misidentfied as muon) are a rather 
small background component. Except for hadronic 
decays, which are negligible, all backgrounds are com- 
mon to the ev-mode, which is suppressed by a factor of 
0(10^^). Thus, the ei^-sample provides a good model 
of the fiv background that has to be corrected only for 
kincmatical and efficiency differences. Including this 



corrected shape in the fit, the total number of fitted 
/ii/-events in data is given by 

A^r; = E<''<'''^'\ (18) 

i=l 

where A'^^*-"'* represents the total number of recon- 
structed signal MC events that were generated in the 
i-th bin of (regardless of the reconstructed n^) 
and uif is the fitted weight of this component. 

The numerical result for Nj^"" is 32100 ±870(stat)± 
1210(syst), that for NJ^l^ is 169 ± 16(stat) ± 8(syst). 
The statistical uncertainties are due to statistics in 
the data signal, the systematic uncertainties due to 
statistics of the data background samples and those 
of the MC samples used. These errors include the 
non-negligble correlations between the ri^ bins. 

As the branching fraction of Ds used for 

the generation of generic MC is known, the branching 
fraction in data can be determined using the following 
formula: 

^rcc 

B{Ds fliy) = GMCcxp ■ ^generated (£»« ^ M^^) , 

(19) 

where 6gencrated(i^. ^ /^i') = 0.0051 and iVGMCcxp ig 

the number of reconstructed ^i^-events in the generic 
MC, weighted according to the fit to data, i.e. 

^GMCexp^^^Z^.^GMC,., (2O) 

where A^^^^'* represents the total number of recon- 
structed /ii^-events filtered from GMC that were gen- 
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FIG. 10: Belle: B{Ds fiu) as a function of n^; final re- 
sult is shown as the green shaded region. For comparison, 
the PDG value and its error is shown as yellow shaded 
region in the background. 



erated in the i-th bin of (regardless of the recon- 
structed n^). 

Figure [TO] shows the branching fraction determined 
in bins of nj^ (using correlated fit results). The re- 
sult is stable within errors in nx ; note that the errors 
shown for the nx bins are the total errors, including 
correlation. The final result is: 

BiDs fiiy) = (6.44 ± 0.76(stat) ± 0.52(syst)) • 10"^ 

(21) 

The statistical error reflects the statistics of the signal 
sample. The systematic error is dominated by statisti- 
cal uncertainties due to background samples from data 
and MC samples (0.29) and the statistical uncertainty 
on 7V^j^^°'^P (0.41). Since the branching fraction is 
determined by calculating a ratio of the signal yield 
to the number of D^-tags, systematics in the recon- 
struction of the tag side cancel; the only remaining 
systematics are due to the tracking and identification 
of the muon, which have been estimated as 2%, con- 
tributing 0.13 to the total systematic error. As a cross- 
check, also the branching fraction for < 6 has been 
determined as (6.54 ± 0.76(stat) ± 0.54(syst)) • 10'^, 
which agrees nicely with the result given above includ- 
ing all available njr bins. 

The decay constant fn^ , using Eqn. [16] and recent 
values from PDG [lB| yields 

fo, = 275 ± 16(stat) ± 12(syst)MeV. (22) 



VI. SUMMARY: OVERVIEW AND 
COMPARISON 

Studies in the charm sector are notoriously difficult 
for experiments running at much higher than thresh- 
old energy; still both Belle and BaBAR present an in- 
teresting variety of results on (semi)leptonic charm de- 
cays. While Belle concentrated on fully reconstructed 
(and consequently tagged) events, BaBAR preferred 
methods with partially reconstructed events, and uses 
tag information only for its Dg IJ-i^ii analysis. 

The advantage of Belle's approach is a very effec- 
tive background suppression and an excellent neutrino 
momentum resolution which can compete with results 
achieved at exp eriments operating at threshold energy 
like Cleo-c [23|. It also allows absolute measurements, 
by normalization to the number of tags. 

However, BaBAR's approach is significantly more ef- 
ficient in terms of event statistics. In the case of the 
D'^ — !■ KeVe analysis, despite higher backgrounds, this 
results in a much better accuracy of measurements. 
For the Dg ijlv^ analysis, the advantage of higher 
statistics is more or less equalled by the disadvantage 
of larger backgrounds, which eventually gives some- 
what larger errors than Belle's. 

In any case, it is a valuable crosscheck to have rather 
different experimental approaches at different experi- 
ments. Table lllll summarizes all results by Belle and 
BABm. 

Within the semileptonic channels discussed in this 
review, only the mode — > Keve has been studied 
by both BABm. and Belle, and can be compared. The 
measured branching fractions agree well within errors, 
larger differences are seen in the form factor measure- 
ments. However, these differences do not exceed 1.3it, 
and could be due to the more complicated systematics 
of the fits involved. The results are much more precise 
than those of previous experiments, and also well com- 
patible with other recent results [H, [2^ . Obviously, 
the untagged, partial reconstruction of BABMi has 
clearly smaller erros, even though it suffers a large un- 
certainty due to the normalizing channel used, which 
is the dominating part of its systematic error. Belle 
does an absolute measurement, but also is clearly lim- 
ited by systematics. Thus in neither case a significant 
further improvement of the measurements can be ex- 
pected with more data accumulated, without also fur- 
ther developing the experimental methods. 

The other mode where results can be compared is 
Ds — > ^Vfji- The results agree very well with each 
other. Here Belle profits from its full reconstruc- 
tion method, and has somewhat, but not dramatically 
smaller errors. Both results are well compatible with 
other recent results, and still compatible with theoret- 
ical predictions, which tend to be lower by « 2.5(Texp, 
but bear some uncertainties as well [30| . In both ex- 
periments, statistical and systematic error are of the 
same order. Considering the fact that part of the 
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TABLE III: Overview of results obtained by BaBar and Belle experiments, compared to theoretical expectations (where 
available); errors are statistical (first) and systematic (second). 



decay mode 



parameter 



D -» Jfez/e BF (%) 

z-expansion, ao 
z-expansion, oi 

z-expansion, 02 
w-poie (GeV/c^) 
a 

/+(0) 

D KtiVi, BF (%) 

jTipoie (GeV/c^) 
a 
/+(0) 

D Trei/e BF (%) 

mpoie (GeV/c^) 
a 
/+(0) 

D Tvniy^ BF (%) 

mpoie (GeV/<?) 

a 

/+(0) 

Da — > (j>eve Ti {mA,v fixed) 
rv {mA,v fixed) 
r2 {mv fixed) 
rv {mv fixed) 
ruA (GeV/c^) (my fixed) 
/ttfp BF (10-^) 
fp, (MeV) 



prediction 



2.112(= moj) 
0.50 ± 0.04 
0.73(3)(7) 

2.112(= moj) 
0.50 ± 0.04 
0.73(3)(7) 

2.010(= mo*) 
0.44 ± 0.04 
0.64(3) (6) 

2.010(= mzj.) 
0.44 ± 0.04 
0.64(3) (6) 



249 ± 3 ± 16 



BaBar result 



3.522 ± 0.027 ±0.079 
2.98 ±0.01 ±0.04 
-2.5 ±0.2 ±0.2 
0.6 ±6.0 ±5.0 
1.884 ±0.012 ±0.015 
0.377 ± 0.023 ± 0.029 
0.727 ± 0.007 ± 0.009 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 
n/a 

0.705 ± 0.056 ± 0.029 
1.636 ± 0.067 ± 0.038 
0.711 ±0.111 ±0.096 
1.633 ± 0.081 ± 0.068 
2.53l°jt±0.54 
6.74 ±0.83 ±0.71 
283 ± 17 ± 16 



BcUo result 



3.45 ±0.10 ±0.19 
n/a 
n/a 

n/a 

1.82 ± 0.04 ± 0.03 
0.52 ± 0.08 ± 0.06 

0.695 ±0.007 ±0.022 
3.45 ±0.10 ±0.21 
included in KeVe 
results shown 

above 

0.279 ±0.027 ±0.016 
1.97 ±0.08 ±0.04 
0.10 ±0.21 ±0.10 
0.624 ± 0.020 ± 0.030 
0.231 ± 0.026 ± 0.019 
included in Treve 
results shown 
above 
n/a 
n/a 
n/a 
n/a 
n/a 

6.44 ±0.76 ±0.52 
275 ± 16 ± 12 



difi'erence 



systematic error is due to the size of control samples 
which will get larger with more statistics, there is some 
room for further improvements once the full data sets 
of the experiments are available. 
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